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Abstract. We present a detailed photoionization model of G29.96-0.02 (hereafter G29.96), one of the brightest 
Galactic Ultra Compact H n (UCHII) regions in the Galaxy. This source has been observed extensively at radio 
and infrared wavelengths. The most recent data include a complete ISO (SWS and LWS) spectrum, which displays 
a remarkable richness in atomic fine-structure lines. The number of observables is twice as much as the number 
available in previous studies. In addition, most atomic species are now observed in two ionization stages. The 
radio and infrared data on G29.96 are best reproduced using a nebular model with two density components: a 
diffuse (n e ~ 680 cm -3 ) extended (~1 pc) component surrounding a compact (~0.1 pc) dense (n e ~ 57000 cm -3 ) 
core. The properties of the ionizing star were derived using state-of-the-art stellar atmosphere models. CoStar 
models yield an effective temperature of ~ 30*i kK whereas more recent non-LTE line blanketed atmospheres 
with stellar winds indicate somewhat higher values, T c s ~ 32-38 kK. This range in T e s is compatible with all 
observational constraints, including near-infrared photometry and bolometric luminosity . The range 33-36 kK is 
also compatible with the spectral type 05-08 determined by Watson fc Hanson (199T| ) when recent downward 



revisions of the effective temperature scale of O stars are taken into account. The age of the ionizing star of 
G29.96 is found to be a few 10 6 yr, much older than the expected lifetime of UCHII regions. Accurate gas phase 
abundances are derived with the most robust results being Ne/S = 7.5 and N/O = 0.43 (1.3 and 3.5 times the 
solar values, respectively). 
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1. Introduction 



The properties of young massive stars are still poorly 



coupled with stellar atmosphere models allows one to con- 
strain the properties of the ionizing stars. 



known i This 



is due m me lad mm massive sfais aie it 
atively'iaie and thai Llieii lifuliiiiu is sliuiL In Lliuii eaily 

stages, massive stars are still embedded in their parental 
molecular clouds and suffer heavy extinction. Our knowl- 
edge of the stellar energy distribution in the visible and 
the ultraviolet of young massive stars therefore relies on 
indirect measurements, such as atomic fine-structure lines 
observed at infrared wavelengths. A comparison of the line 
fluxes with predictions of detailed photo-ionizing codes 



Recently, Peeters et aL (2002) and Martin-Hernandez 
ct al. (2002| ), hereafter Papers I & II, presented the results 
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of an infrared spectral survey of 34 galactic compact H n 
regions based on complete ISO grating spectra. For most 
of the sources, the ISO data display a remarkable richness 
in spectral lines including all the fine-structure lines of N, 
O, ne, S, and Ar. Using additional data from the literature, 
the elemental abundances and their variation across the 
Galactic disc were derived (Paper II). 

Some of the H n regions in the sample are well known 
and were studied in detail at other wavelengths. Those ad- 
ditional data provide useful contraints to derive the prop- 
erties of the ionizing star and to fine tune the elemental 
abundance estimates. One such source is the compact H n 
region IRAS 18434-0242 (G29.96-0.02, herafter G29.96). 
This compact source has one of the richest ISO spec- 
trum of the entire sample (Paper I) and is one of the few 
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compact H n regions for which the ionizing star has been 
identified and characterized through direct infrared spec- 
troscopy (Watson & Hanson 1997; Hanson ct al. 2002). 

The availability of high-quality data on G29.96 to- 
gether with the recent developments of the models of stel- 
lar atmosphere of massive stars prompted us to make a 
detailed photionization model of this source with the aim 
to further constrain the nature of the ionizing stars and 
to derive accurate elemental abundances. The paper is or- 
ganized as follows: the observational properties are sum- 
marized in Sect. 2; Sect. 3 describes the photoionization 
code, the input stellar parameters and the methodology; 
Sect. 4 presents the results of the best model which are 
discussed in Sect. 5; the conclusions are given in Sect. 6. 

2. G29.96-0.02: Observational facts 



The compact Hn region G29.96 is one of the brightest ra- 
dio and infrared sources in our Galaxy. Its morphology is 
a classical example of a com etary-like compact H n region 
( Wood fc Churchwell 1989b) in inter action with a molec- 
ular cloud (e.g., |Pratap et al. 199S| , hereafter PMB99). 
G29.96 has been studied in detail at infrared and radio 
wavelengths and the following summarizes the main re- 
sults together with the essential properties of this compact 
H II region. 



2.1. Observations of fine-structure lines 

The fine-structure lines observed by the SWS and LWS 
spectrometers are tabulated in Paper I where their ob- 
served fluxes associated error bars are given as well as the 
entire ISO spectrum of this source. In addition, 11 hydro- 
gen recombination lines were detected in G29.96 and were 
used to determine the interstellar extinction (Paper II). 
We will use Bra in the modeling process (see Sect, hi 
and Table [l]) as this is the brightest H I line and one of 
the least affected by extinction. The ISO lines corrected 
from the interstellar extinction (using Ak= 1-6 and the 
extinction law derived in Paper II) are given in Column 2 
of Table |. 

Infrared fine-structure lines have previously been ob- 
served with KAO by Hcrtcr et al. (198l|) , who measured 
the lines [Ar il] 6.98/mi, [Ar in] 8.98^m, [S in] 18.7/zm and 
[Nell] 12.8/im, in agreement with the ISO values within 
20 %. 

Similarly, (Bimpson ct al. (1995| ) reported obser- 
vations of [S in] 33.6/im, [Nem] 36.0/im, [O ill] 51.8/mi, 
[Nm] 57.3/zm and [O ill] 88.3/xm differing from the ISO 
values by less than 10 %, except the [Siii]33.6^m and 
the [O in] 88.3/^m fluxes which are both ~ 30 % higher in 
our sample. The difference of aperture sizes and pointing 
between KAO and ISO can partly be responsible of these 
discrepancies. 



M aps o f [Neil] 12.8/mi were obtained by Lacy et al 
1982]) and [Watarai et al. (1998| ) . These authors found an 



2.2. Radio observations: core/halo structure and 
constraints on the He ionization structure 

G29.96 has been observed at radio frequencies using var- 
ious spatial resolutions. At 2 cm the observed flux densi- 
ties range from 2.7 to 4.6 Jy (|Wood fc Churchwell 1989b 



Afflcrbach et al. 1994| : [Fcy ct al. 1995|; |Watson et al. 1997|) 
at 6 cm from 1.4 to 3.6 Jy (|Wood fc Churchwell 1989b 



Afflcrbach ct al. 1994), and at 21 cm from 0.9 to 2.6 Jy 



flClaussen fe Hofner 1995| ; [Kim fc Koo 200l| ). As usually 
the case for UCHII regions, the source presents a diffuse 
emission. Therefore the highest spatial resolution obser- 
vations may miss part of the radio flux density. We have 
adopted the following values: 3.9, 3.4 and 2.6 Jy at 2, 6, 
and 21 cm respectively, favoring the highest values to take 
into account the diffuse emission. The resulting number of 
Lyman continuum photons (1.8 x 10 49 , cf. below ) is higher 
than the value of 4 x 10 48 derived in Paper II, but this last 
value was obtained for a uniform electron density of 10 4 
cm~ 3 and a size of 7 arcsec. 

G29.96 is characterized by a strong edge-brightened 
core/" head", with a low surface-brightness "tail" of emis- 
sion trailing off opposite the bright edge. Due to this ex- 
tended emission , part of G29.96 cannot be s trictly called 
Ultra Compact. Wood fc Churchwell (1989b| ) obtained n e 
= 8.5xl0 4 cm~ 3 in the arc and a n e 5-10 times lower 
in the tail of the nebula, while Afflerbach et al. (1994) 



estimate n e = 5.x 10 cm 3 in the leading arc and n e 



= 2.xl0 4 cm 3 in the tail. |Simpson et al. (1995 ) ob- 



tained n c = 1500 cm~ 3 from the [O in] 51.8/ 88.3^m 
lines ratio. From the set of the ISO observables avail- 
able for G29.96, 3 line ratios can be used as density 
diagnostics: [O in] 51.8/ 88.3^m, [S in] 18.7/ 33.6^m and 
[Ne in] 15.5/ 36.0/^m. Unfortunately, the two last ratios 
suffer large calibration uncertainties (25 % at 1 <r error, 
see Paper I) and therefore only [Oni] 51.8/88.3/im (here- 
after r[om]) can be safely used to derive the gas density. 
From rrom] = 2.4 an electron density n e ~ 800 cm~ 3 is 
derived (Paper II). 

This apparent discrepancy between the densities deter- 
mined from fine-structure lines of oxygen and from radio 
observations have already been pointed out by Afflerbach 
[et al. (1997 ) who suggested a core/halo description of 
G29.96. Faint diffuse halos are commonly obs erved in the 
radio continuum maps of UCHII regions (e.g., Garay et al 



1993 ; [Fey et al. 1995|; |Afflerbach et al. 1996| ; [Kurtz et al 
1999[ ). Recently, [Kim fc Koo (200l| ) found extended cmis- 



integrated flux of 84 % and 64 % of our value, over a size 
of 10" x 10" and a diameter of 30" respectively. 



sion at 21 cm linked to the bright spot of G29.96. Bra im- 
ages ( Lumsdcn fc Hoare 1996| ; |Watson et al. 1997 ; Watson 
|fc Hanson 1997| , PMB99) also support this morphology. 

The observational evidence of a core/halo morphology 
has led us to model G29.96 with two components as out- 
lined in Sect. |3.3| . 

Kim & Koo (2001) have observed the He76a and H 
radio recombination lines in G29.96. For various positions, 
including the region considered here, they obtain a He + 
abundance of ~ 0.068-0.076. Assuming a normal helium 
abundance of 0.1 this implies that He is predominantly 
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singly ionized helium in G29.96. This provides a useful 
constraint on the temperature of the ionizing source (see 
Sect. O). 



3.1. Photoionization code 

The models are performed using the detailed photoioniza 
tion code NEBU (Morisset fc Pcquignot 1996; Pcquignot 



2.3. Heliocentric and galactocentric distances 

With an average global radio recombination lines LSR 



ct al. 2002 ) which consistently computes the line fluxes 



velocity of 95 km/s (Afflerbach et al. 1994), adopting a 
galactocentric radius of the Sun of 8.5 kpc, a kinematic 
heliocentric distance for G29.96 between 5.5 and 9.5 kpc 
is derived using a standard galactic rotation curve with a 
rotation speed of 220 km/s at the Suns position. Previous 
investigators assumed the average between the near and 
far distances. 

The extinction along the line of sight to G29.96 was 
studied by PMB99 using galactic Hi and CO surveys. 
They found Ak ~1 at a distance of 5 kpc, and Ak ~3 at a 
position corresponding to the tangent point (at about 7.5 
kpc). These findings provide a crucial argument to con- 
sider that the near distance should be more appropriate. 
Hereafter, we will adopt an heliocentric distance of 6 kpc 
(+1.0,-0.5) for G29.96Q, implying R G =4.5 kpc (+/-0.3). 

2.4. Ionizing star(s) 

A detailed search for stars embedded in the H II region and 
the adjacent molecular hot core has been performed in the 



near-infrared (Lumsdcn & Hoarc 1996; Watson et al. 1997|; 



Watson fc Hanson 1997j , PMB99). |Watson et al. (199 
and PMB99 have in particular revealed the existence of 
a cluster of about 18 OB-type stars, or their progenitors, 
embedded in the cloud. The same authors have convinc- 
ingly identified the bright star at the center of the arc of 
radio emission as the exciting star, or at least as the pri- 
mary source of ionization. There is no evidence for an in- 
frared excess in the K band, suggesting that any remaining 
disk should be optically thin, and therefore that the star 
is no longer accreti ng mass. From near infr ared spectro- 
scopic observations, Watson fc Hanson (1997 ) were able to 
constrain the spectral type of the ionizing source. A spec- 
tral type between 05 and 08 (and no constraint on the 
luminosity class) w as found. In contrast, a recent study by 
Kaper et al. (2002 ) reports a spectral type as e arly as 03 



From their spectroscopic and photometric data Watson fc 



Hanson (1997 ) and Watson et al. (1997 ) already deduced 
an evolutionary age of about 1-2 x 10 6 years for a single or 
binary star, in apparent contradiction with the estimated 
age of the UCHII region (~ 10 5 yr). 



The overall SED of G29.96 derived by Watson et al 



(1997 ), the near-IR photometry, and spectral types pro- 



vide important constraints on the fundamental properties 
(T ff, luminosity, age) of the ionizing source. A detailed 
discussion is given in Sect. 5.7. 



3. Photoionization modeling 



1 At 6.0 kpc, 1" = 0.029 pc = 8.95xl0 16 cm 



without any hypothesis on the ionization structure of the 
gas, especially without using ionization correction factors 
(zc/s). However, it does require assumptions about the 
geometry, density, and pressure structure of the nebula. 

The computation is performed in a spherical geometry 
and at each radius from the central ionizing source, the 
electron temperature, the electron and ions densities, and 
the line emissivities are determined solving the ionic and 
thermal equilibrium equations. The inputs for the model 
are the description of the ionizing flux, using e.g. an effec- 



tive temperature and a luminosity (see Sect. 3.2), and the 
gas distribution (assuming e.g. constant pressure through 
the shell) with a set of abundances. 

The elements taken into account and for which lines 
fluxes are predicted are: H, He, C, N, O, Ne, Mg, Si, S, 
CI, Ar, Fe, Ca and Ni. Self-absorption effects of the radio 
continuum are computed in a spherical geometry approx- 
imation. 

Absorption of incoming and diffuse photons by dust is 
considered, the number ratio of dust grains over the num- 
ber of hydrogen atoms being a parameter of the model. 
The optical properties of dust grains were made available 
to us by Ryszard Szczerba (private communication). The 
adopted dust composition is 50 % " astronomical" silicates 



and 50 % graphite (Draine & Laor 1993). No quantum 
heating by very small grains are taken into account in the 
version of NEBU used for this work. 

3.2. Atmosphere models 

The ionizing photons distribution is taken from the re- 
cent CoStar atmosphere models, which in clude non-LTE 
effects, line blanketing, and stellar winds ( Schaerer fc de 
Kotcr 1997). The CoStar atmospheres have been com- 



pared to earlier predictions (Kurucz models and plane par- 
allel non-LTE atmosphere models). Their implications on 
nebular studies have also been discussed by Stasinska & 
[Schaerer (1997 ) and Schaerer fc Stasinska (1999] ). 

The predicted ionizing fluxes have been "tested" by 
the following indirect means: 

1) the CoStar SEDs reproduce well the ionic Ne++/0 ++ 
abundance ratios in Galactic Hn regions observed 
from optical and I R data, thereby solving the s o-called 
[Ne in] problem (|Stasihska fc Schaerer 1997^ . This 
supports the important increase of the ionizing spec- 
tra predicted at energies 2; 35-40 eV. 



The atmosphere models calculations of Sellmaier et al 



(1996) have not been confirmed by the Munich group. The 



spectra from their latest grid of atmospher e models includ- 



ing n on-LTE line blocking and blanketing (Pauldrach et al 



2001 ) predict [Ne ill]/ [Ne n] ratios in r ough, though somew hat 
marg inal agreement with observations ( Gtiveon et al. (2002 ) , cf . 
Sect. 
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Fig. 1. Emission lines maps of the best model, in grey linear scale. The two componets are presented with angular 
size related to their corresponding covering factors. Contours for the ISO SWS beam profile are superposed, for 
transmissions of 20, 40, 60, 80 % of the flux. In the four latest maps, the LWS aperture size is of order of the images, 
and not shown. Note that the SWS aperture is not always centered. 



2) the observed increase of the ratio of the He over 
H° ionizing photons between stellar effective temper- 
atures of ~ 30000 and 40000 K is well reproduced 
( |5chacrer 2000 ). This quantity, sensitive to the ion- 
izing flux above 13.6 and 24.6 eV, is constrained by 
optical He I recombination line measurements in the 
sample of Kcnnicutt ct al. (200C| ) of Hll regions with 
known stellar content. 

3) Tailored photoionization models for two nebulae with 
exciting stars of spectral types 03-04 and 07 respec- 
tively have been constructed by pey et al. (200"0 ) using 
CoStar spectra. While overall the spectra are found 
to yield a good agreement with the observed nebu- 
lar line ratios for all objects (confirming also the first 
point), there is some indication of an overprediction 
of the spectrum above ~ 35-40 eV for the early-type 
object (DEM L323). We note, however, that for DEM 
L243, this discrepancy is not found if a cooler spec- 
trum is adopted (approximately model B2 instead of 
C2 preferred by Oey et al.), as would be expected for 
the region whose ionizing flux is likely dominated by 
the Q7If supergia nt instead of the 07V((f)) dwarf (cf. 
Vacca et al. 1996|). 



For the range of stellar temperatures of interest here 
(T e g ~ 28-34 kK; cf. below) only 1 and 2 are directly 
applicable here. Indeed, at these relatively cool effective 
temperatures the uncertainties in the CoStar models are 



expected to be the largest, as stressed by Schaerer & de 



Koter 1 1997). Although the constraints on the atmosphere 
models are still rather limited, we conclude from the above 
comparisons that the adopted model atmospheres should 
provide a reasonable description of the ionizing fluxes. 
Comparisons with other model atmospheres are discussed 
in Sect. 5.6. 

From the 27 CoStar models available on the Web 



distributions was constructed. For any effective tempera- 
ture T e g and luminosity L in the range covered by the 
CoStar models, an interpolation is performed between the 
four nearest models of this grid, using the square inverse 
of the distance in the log T e ff-log L plane to determine 
the weights of the four CoStar models. The four CoStar 
models are first divided by the corresponding blackbody 
spectra, then averaged using the weights previously deter- 
mined and the result is finally multiplied by the blackbody 
spectrum of the desired T g and L, avoiding the diktat of 
the most intense CoStar model. 



3.3. Two-component model 

In order to reproduce both rrom] an d the ra dio fl ux den- 



( (Bchacrcr fc dc Koter 1997 ), a finer grid of spectral energy 



sities which imply a higher density (see Sect. |2.2| ), a two- 
component model is used to describe G29.96. The two 
components differ by their densities; the lower and higher 
density components are named component 1 and 2 respec- 
tively. 

The present model consists of a simple linear combina- 
tion of two independent runs of NEBU in a spherical case 
and under isobar approximation. Both components are as- 
sumed to be radiation bounded. The gas is supposed to 
be homogeneously distributed in each component, with an 
inner radius to reproduce an empty cavity The coefficient 
applied to each component to obtain the fluxes of the emis- 
sion lines is the covering factor, i.e. the angular size over 
4 7r of each component as seen by the central source. The 
sum of the two covering factors is set to one, i.e. we do not 
consider that photons could escape from the Hll region. 
In this model, there is no diffuse radiation exchanged be- 
tween component 1 and 2, and no effects of shielding of 
one component by the other is taken into account as: com- 
ponent 2 has a very small geometrical thickness compared 
to component 1, and component 1 is very optically thin at 
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radio wavelengths. A two-component model must be seen 
as an approximation, describing the two first moments of 
a certainly more complicated gas distribution. 

The most important effect of this 2-density medium 
is in the localization of the lines emission, which depends 
on their critical densities. For example, the [O III] and 
[N III] lines are emitted only by the lower density part of 
G29.96, because these lines are collisionally de-excited in 
the densest region (see Paper II for the critical densities 
of all the lines). 

3.4. The ISO beam profiles 

Component 1 of the nebula will have a larger geometric 
extension than component 2 due to its lower density. The 
beam sizes of the ISO instruments, which vary from 14- 
33 arcsec for the SWS to ~ 80 arcsec for the LWS, have 
a drastic impact when comparing/di viding different cmis - 
sion lines. The beam profiles, from Garcia-Lario (1999| ), 



1 0.0 



are used to apply a correction to the predicted line fluxes. 
For each component and line, an intensity map is gener- 
ated by projection of the line emissivity on a sky plane. 
An "ISO" mask is then applied, according to the detector 
beam corresponding to the line wavelength. 

Fig. ^ shows the intensity maps obtained for our best 
model (presented in Sect. |J) for 8 atomic lines. Contours of 
the corresponding ISO SWS apertures are superimposed 
on the image. For the four last images corresponding to 
lines observed by LWS, the aperture size is larger than 
the image. For most of the lines, only the core is visible 
with the adopted linear intensity scale. Note that, despite 
the absence of visible extended emission, for most of the 
lines the contribution of the extended part is about that of 
the core. The low density extended component is clearly 
seen for the low critical density lines for which the contri- 
bution of the core is very low (the [Nil], [NIII], [OIII] 
and [S III] lines). The effect of neglecting the finite aper- 
ture size is clearly illustrated. Note also that the profil e of 
the S WS aperture, as described by Garcia-Lario (1999). is 



1 .0 



0.1 



26 




30 32 
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not sy mmetrical. We did not attempt to exactly adjust the 
orientation of these profiles to the observations of G29.96, 
as the effect of the asymmetry in the profiles is negligible 
compared to the effect of departure from spherical sym- 
metry of the object itself. 

4. Results 

4.1. Convergence procedure 

Table |] lists the parameters of the best two-component 
model and the observable used to constrain each param- 
eter. In the following we describe the convergence proce- 
dure of this model. 

The [O in] lines are collisionaly de-excited in compo- 
nent 2 and trace mainly component 1, while the other 
density diagnostic lines are emitted by both components 
of the nebula. Therefore the [O in] 51.8/ 88.3/im ratio is 
used to constrain the density of component 1. The den- 



Fig. 2. Variation of the four [X ,; ]/[X* +1 ] ratios with the 
effective temperature. The ratios correspond to: Ar (solid) 
, S (dot), Ne (dashed) and N (dot-dashed) and have been 
divided by their corresponding observed values. 

sity of component 2 is constrained by the 6 cm radio flux 
density. The ratio of the two covering factors is determined 
by fitting the Bra line flux, the sum beeing fixed to 1. 

Fig. | shows the four available [X 1 }/[X l+1 ] line 
ratios (divided by the corresponding observed val- 
ues) versus the effective temperature of the ion- 
izing star, namely [Aril] 6.98/im/[Anii] 8.98/im, 
[S m] 18.7/im/[Slv] 10.5/xm, [Neil] 12.8^m/[Neni] 15.5/xm 
and [Nil] 121.7/mi/[Nm] 57.3/xm. The number of ionizing 
photons is kept constant by adjusting the number of 
ionizing stars with fixed luminosity, while the effective 
temperature is varied. The models used for this plot are 
all 2-component models. In the range of T c g considered 
here (27 to 35 kK) the sensitivity of those ratios to the 
effective temp erature is very high (see also Stasihska & 
Schaerer 1997[ and Paper II). Within a range of ~ 2500 K, 



the four ratios provide the same constraint on the effective 
temperature. However, the [Am] 6.98/xm/[Arni] 8.98/im 
ratio gives a higher effective temperature than the three 
other diagnostics (see Sect. 5.6). We will therefore use 
the [Nil] 121.7/xm/[Nlll] 57.3/xm ratio to determine the 
effective temperature, as it is: 1) independent of the ISO 
beam size since both lines were observed with LWS, and 
2) emitted only by component 1, since it is collisionaly 
de-excited in component 2. An effective temperature of 
~ 30l^ kK is then derived. 

The total luminosity is the product of the CoStar 
model luminosity by the number of stars, and the num- 
ber of ionizing photons is adjusted to reproduce the radio 
flux density at 2 cm. A range of stars with different spec- 
tral types is presently not considered. With the available 
CoStar models, there is an upper limit for the star lu- 
minosity beyond which no more models are available and 
which corresponds to the post main sequence or Wolf- 



6 



C. Morisset et al.: A photoionization model of G29.96-0.02 



Table 1. Parameters of the best model and derived gas properties. For some parameters, the observable used for the 
main constraint is given. 

Parameters common for both components 



Effective temp (kK) 

Luminosity in log(L/Lo) 

Number of stars 

Inner radius (10 17 cm) 

Dust/gas ratio 

[He]/[H] (in numbers) 

[C]/[H] 

[N]/[H] 

[0]/[H] 

[Ne]/[H] 

[S]/[H] 

[Ar]/[H] 



29.7 
5.86 1 
1.5 
3. 

10" 5 
0.1 

1.00xl0~ 4 
1.97xl0 -4 
4.55xl0~ 4 
1.70xl0~ 4 
2.25xl0" 5 
5.00xl0~ 6 



from [N n] 121.7^tm/[N ill] 57.3/zm 

from borderline in CoStar models grid (see Sect. 5.8) 

from 2 cm flux density (see Sect. 5.5) 

from imaging 

from ISO continuum emission 



from [Nil] 121. 8^m + [Nm] 57.3/um 
from [O in] 51.8 + 88.3/jm 
from [Neil] 12.8/^m + [Nem] 15.5/wn 
from [S ill] 18.7Atm + [S iv] 10.5^m 



Parameters for: 


Component 1 


Component 2 


Inner \\h (cm -3 ) 


640. from rjom] 


52000. from 6 cm flux density 


Covering factor 


36 % and 


64 % from Bra 


Filling factor 


1.00 


1.00 


Properties for: 


Component 1 


Component 2 


Mean n e (cm~ d ) 


680. 


57000. 


Constant pressure (cgs) 


9.7xl0~ 10 


l.lxlO -7 


Thickness (10 17 cm) 


27.2 


0.15 


Mass (M ) 


29. 


0.73 


Mean T H n (K) 


5520 


7230 


inner U 


8.8xl0 _1 


5.7xl0" 3 



This luminosity multiplied by the number of stars leads to 1.8 x 10 49 ionizing photons (NL yc 



Rayet star regime. The degeneracy of the CoStar m odel 
luminosity by the number of stars is discussed in Sect. 5.5 . 

Once the above set of parameters is fixed, the abun- 
dance of each species is derived by fitting their lines fluxes. 
As most of the parameters have feedback effects, an iter- 
ative process must be applied. 

To summarize, we have 17 free parameters (listed in 
Table [l]), namely: effective temperature and luminosity of 
one star, number of stars, densities of both components, 
ratio of the covering factors, inner radius of the empty cav- 
ity, dust to gas ratio, filling factors of both components, 
and 7 abundances. On the other hand, we have 16 observ- 
ables plus the morphology of the source as given by the 
radio maps. 

Some parameters (the inner radius, the dust/gas ratio, 
the filling factors) cannot be precisely constrained and are 
set to a reasonable value. Change in their values are dis- 
cussed in Sect. |[ The abundance of helium and carbon 
does not have any effect on the results of the model, if 
remaining within reasonable values. 

We finally stay with 11 free parameters to be con- 
strained by 17 observables. In other words, there are 
6 observables that are not used to build the model 
and are therefore entirely predicted, namely: the three 
[X J ]/[X 4+1 ] ratios for Ne, Ar and S, the two [S hi] 33.6/^1 
and [Nem] 36.0^m line fluxes (not used because of their 
large calibration error, see Paper I and Paper II), and the 
21 cm continuum flux density. 



4.2. Results and predictions 

Table || lists the observations of the infrared emission lines 
and the radio continuum flux densities together with the 
results of the best model. The contributions of the two 
density components are given separately. Most of the pre- 
dicted lines fluxes and radio flux densities agree with the 
observations to within the uncertainties of the measure- 
ments. The few predictions which are off by a larger fac- 
tor are the results of well understood factors which will be 
explained hereafter. 

The three [X*]/[X l+1 ] ratios not used to determine 
the effective temperature agree within ~ 1 kK with 
the [Nil] 121.7/xm/[Nm] 57.3/im ratio (see Fig. |). The 
[SIV] 10.5/^m line falls in the absorption feature due to 
silicate. This has been taken into account when using the 
attenuation law described in Paper II and the model pre- 
diction is in very good agreement with the observation. 
For [Nelll] 15.5/xm, we can suspect an overprediction of 
the line fl ux due to the use d of the CoStar models, as dis- 
cussed by |Oey et al. (2000 ) and in Sect. [T^ (see also Sect. 



Note that the predicted 21 cm flux densit y is lower 
than the value observed by Kim fc Koo (2001 ), who re- 
ported complex extended radio emission toward G29.96. 
Part of this diffuse emission could be due to gas ionized 
by members of the stellar cluster other than the main ion- 
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Table 2. Observations and model predictions. The model fluxes for each component arc already multiplied by the 
ISO SWS aperture correction factor (given in parenthesis for component 1, component 2 being small enough to be 
entirely seen by both ISO spectrometers) and by the covering factor corresponding of each component. 



Line 




Line fluxes (10 -18 W/cm 2 ) 




Model/Observation 




Observations 


1 Model 


Component 1 


Component 2 




Hl4.05/im 


11.4 ±0.4 


11.4 


.900 (.16) 


10.5 


1.00 


[Arii]6.98/im 


27.1 ±3.0 


36.3 


1.34 (.08) 


34.9 


1.34 


T A r TTTl & Oa/jTvi 


9fl ^ -4-1 A 


13.1 


2.21 (.23) 


in q 


U.D4 


[o IIIJ lo. ( /im 


a A *7 _1_1 o 

04./ ±1.2 


55.5 


25.0 (.21) 


on a 
3U.4 


U.86 


re tttI qq (?, 
[o IIIJ oo.D/im 


/I A q2 _i_q n 

44. o ±O.U 


23.8 


20.8 (.18) 


z.9o 


n K a 2 
U.04 


S Tvl 10 5//m 


7 90 ±0 8 


8.95 


7.70 (.51) 


1.25 


1.13 


[Neil] 12.8/tm 


127.5 ±9.0 


97.8 


6.08 (.13) 


91.7 


0.77 


[Ne ill] 15.5/mi 


34.0 ±1.3 


42.5 


15.9 (.43) 


26.6 


1.25 


[Nein] 36.0/im 


3.28 2 ±1.0 


2.37 


1.18 (.38) 


1.18 


0.72 2 


[Nil] 121.8/im 


4.66 ±0.2 


4.65 


4.50 


.150 


1.00 


[N ill] 57.3pm 


23.0 ±0.7 


23.1 


22.6 


.540 


1.00 


[0 in] 51.8/Ltm 


47.4 ±0.2 


47.6 


45.6 


2.08 


1.00 


[O in] 88.3/im 


19.7 ±0.4 


19.6 


19.4 


.210 


0.99 


Wavelength (freq.) 




Continuum flux density (Jy) 






2 cm (15 GHz) 


3.90 ±0.5 


3.90 


1.34 


2.57 


1.00 


6 cm (5 GHz) 


3.40 ±0.2 


3.38 


1.53 


1.84 


0.99 


21 cm (1.4 GHz) 


2.60 ±0.2 


1.88 


1.60 


.282 


0.72 



1 The line fluxes (from Paper I) are corrected for the interstellar extinction using Ak = 1.6 and the law described in Paper II. 

2 Line measured with the SWS band 4 detector, see Paper I. 



izing star of G29.96. This low excitation gas is not taken 
into account in the present model. 

5. Discussion 

The detailed photoionization model of G29.96 reproduces 
with good accuracy most of the atomic fine-structure line 
fluxes and radio flux densities. It allows one to derive the 
elemental abundances in the gas phase and the properties 
of the ionizing star(s). In the following, we will investi- 
gate how some of the less constrained parameters influ- 
ence the results and discuss the reliability of the derived 
abundances. 

5.1. EfFect of the uncertainties of the observed line 
fluxes on the model parameters 

The observed line fluxes are known with 10 to 20 % ac- 
curacy (Paper I). The effect of these uncertainties on the 
model are not always linear. For example, changing rro in i 
by ±20 % changes the electron density of component 1 
by On the other hand, as the stellar effective tem- 

perature diagnostics are extremely sensitive (as shown in 
Fig. ||) , any change by ±20 % in any of these diagnostic 
line fluxes will have virtually no effect on the determina- 
tion of the stellar temperature. 

Concerning the number of ionizing photons, the prod- 
uct of the number of stars by the stellar luminosity is 
directly proportional to the 2 cm flux density. 

One of the most critical parameters is the contribution 
of each component to the total line fluxes. Once the den- 



sity of component 2 is determined from the 6 cm flux den- 
sity, the ratio of the covering factors is derived by fitting 
the Bra line flux. However, this line is sensitive to atten- 
uation and aperture size corrections. As given in Table ^[ 
the nitrogen and oxygen lines are emitted mainly (96 to 
99 %) by the diffuse component where only 10 % of the 
hydrogen lines emission is observed. Decreasing the ob- 
served line flux of Bra by 10 % increases the contribution 
of component 1 from 36 % to 48 % and the density of com- 
ponent 2 from 5.2 xlO 4 to 9.0xl0 4 cm~ 3 . The abundances 
of N, O, Ne, S, and Ar change by -25, -30, ±10, 0, and 
±20 % respectively. 

5.2. Filling factor and components geometry 

The filling factor allows one to artificially increase the geo- 
metrical thickness of the ionized gas. The geometry affects 
both the low and high ionized species if the thickness of 
the nebula is of the order of its radius. 

As component 1 represents the diffuse g filling 
factor of 1.0 seems appropriate (the predicted extension 
of the gas is 35", i.e. comparable to the size of the 
observed surrounding molecular gas). Lowering this fill- 
ing factor to 0.5 has an effect on the lines mainly pro- 
duced by component 1, i.e. [Nil] 121. 8^m, [NIII] 57.3/im, 
[OIII] 51.8, 88.3/im, [Sm] 33.6/xm, and [SIV] 10.5/im. The 
ratio r[oin] remains the same while the [N II] 121.8/im/ 
[NIII] 57.3/xm ratio increases by about 15 %. A small in- 
crease of the effective temperature from 29.7 to 30.1 kK is 
enough to recover the observed ratio. After the whole con- 
vergence process is performed, an increase of the N and 
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O abundances of about 15 % is found. Furthermore, as 
the geometrical thickness of component 1 increases up to 
45" , the effects of the finite size of the ISO SWS beam are 
amplified (the [SIV] 10.5/im line significantly decreases). 
No value lower than 0.5, implying greater geometrical ex- 
tension, would be acceptable. 

For component 2, changing the filling factor from 1.0 
to 0.5 increases the thickness by a factor of about two. 
No obvious effect is found on the line fluxes, but the radio 
flux densities are affected because the self absorption is de- 
creasing with the filling factor. The predicted 6 cm value is 
then higher than the observed value by 8 %, and we have 
to change the hydrogen inner density of component 2 from 
5.2 to 9.2 xlO 4 cm~ 3 to recover it. The geometrical thick- 
ness of component 2 finally decreases from 1.5 to l.OxlO 16 
cm, the changes due to the filling factor being approxima- 
tively compensated by the increase of density imposed by 
the 6 cm flux density constraint. However, the line fluxes 
and the element abundances do not change significantly. 

Decreasing further the filling factor of component 2 to 
a value of 0.1 leads to a different behavior. The hydrogen 
density needs to be increased to 5.5xl0 5 cm -3 in order 
to recover the optical thicknesses of the radio continuum 
at various frequencies. Such a high density implies a col- 
lisional de-excitation of some lines in component 2 (see 
Paper II for the critical densities of all the lines). Finally, 
the abundances are greater than those given in Table | 
by 4, 7, 77, 83, 111 % for N, O, Ne, S, Ar, respectively. 
Oxygen and nitrogen are not very affected as these lines 
are emitted in component 1. The geometrical thickness of 
component 2 becomes 1.6xl0 15 cm. We could interpret 
this model as a distribution of very dense, small clumps 
embedded in the low density medium. 

In summary, modifying the filling factor leads to a 
change of the geometry of the ionized gas. The main ef- 
fect is on the self- absorption of the radio free- free emission; 
changing the filling factor is the same as changing the opti- 
cal depth at the different radio frequencies, in other words, 
changing the proportion of the gas seen tangentially with 
respect to the amount of gas seen radially. 

Afflerbach et al. (1994) derived a filling factor between 
0.03 and 0.4 combining the emission measure obtained 
from the continuum flux density with the local density ob- 
tained from the line-to-continuum ratio. They found high 
values for the density (some 10 4 cm~ 3 ) and they consid- 
ered the gas as included in a sphere: "assuming that the 
line-of-sight depth is equal to the angular diameter from 
the continuum images" . In our case, the dense gas is lo- 
cated in a shell (in which the filling factor is ~ 1.0) with a 
thickness 1/20 of its radius, leading to a total filling factor 
for the sphere of 0.13, compatible with the value obtained 



by Afflerbach et al. (1994). In the model presented in this 
paper, the gas is distributed in a shell, at fixed radius from 
the ionizing source. A more complex model could be con- 
structed with a distribution of clouds at various radii, but 
the new free parameters introduced in such a model could 
not be constrained by any available observable. 



5.3. Role of the inner radius 

We fixed the inner radius of the Hn region to 3 x 10 17 
cm, correspond ing to 3", about the radio core size (e.g., 
Fey et al. 199£ ) . This is also virtually the outer radius of 
the dense component, as the geometrical thickness is 1/20 
times the inner radius. Lowering the value to e.g. 10 17 
cm will require to decrease the density of component 2 to 
2.2 x 10 4 cm -3 in order to recover the radio flux densities 
break between 2 and 6 cm. The geometrical thickness of 
component 2 is now ~ 2 x 10 17 cm, still compatible with 
imaging observations. Oxygen and nitrogen abundances 
must be decreased by 10 %. As the density of compo- 
nent 2 decreases, some lines previously de-excited by col- 
lisions in component 2 are now emitted: [Nem] 15.5/im 
and [S iv] 10.5/xm are predicted to be 3 and 5 times higher 
than the observed values. 

Increasing the inner radius for component 2 will lead 
to an increase of the density of this component to recover 
the radio break. As the geometrical thickness will then de- 
crease to less than one percent of the radius (it was 5 % 
for the adopted model, see Table [l]), the self absorption 
at 6 cm does not increase anymore. This is mainly due 
to the spherical geometry we used for the model. A more 
complicated geometry than such a thin shell could lead to 
more self-absorption. We then relax the constraint of the 
6 cm flux density and perform a model where the 6 cm 
predicted flux density can be higher than the observation 
by some 10 %. With a density of 8 x 10 4 cm~ 3 for compo- 
nent 2, and an increase or N, O, Ne, S, and Ar abundances 
of 12, 12, 24, 30, and 10 % respectively, the result is close 
to the adopted model. 



5.4. Role of dust 

The effect of adding dust in the H n region is to increase 
the absorption of ionizing photons and to change the shape 
of the "apparent" ionizing spectrum. We compared the 
emitted spectra of the dust to the observed infrared con- 
tinuum, at wavelengths lower than 20 /xm. At longer wave- 
lengths, the emission is dominated by cold dust from the 
PDR and the outer H I region, which is not modeled by the 
photoionization code. We check that, whatever the dust 
type used (i.e. "astronomical" silicates, olivines, amor- 
phous carbon or graphite), the modeled emission does 
not exceed the observational data for any concentration 
of dust lower than 5 x 10 -4 (in mass, relative to hydro- 
gen). This represents an upper limit of the amount of dust, 
as part of the 5-20 /im emission can be due to high tem- 
perature PAH's present in the PDR and behind. With this 
amount of dust, the models have to be adjust by changing 
the number of stars from 1.5 to 1.6 and the stellar effec- 
tive temperature from 29.7 to 29.4 kK, without changes 
in the abundances. In this "maximum" dust model, 8% of 
the incoming energy is absorbed by dust in the H n region, 
compared to 11% by ions. The relatively small amount of 
dust derived b y our modeling is in agreem ent with the re- 
cent results of Aannestad fc Emery (2001 ) who found that 
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Fig. 3. HR diagram showing the locations of the available 
CoStar models (crosses) and the four models discussed in 
text (filled diamonds with corresponding number of stars). 
The dotted lines show solar metalicity isochrones for ages 
0, 1.6, 2. 8, 3.5 and 4.0 Myr from left to right, from the 
tracks of [Meynet et al. (199^ ). 



dust in the ionized region of S 125 is severely depleted. 
A more detailed study of the dust emission observed in 
G29.96, including the Hi region, is postponed to a future 
paper. 

5.5. Number of stars versus luminosity and age 

The number of ionizing photons is constrained by the 2 cm 
radio flux density. It is a degenerated parameter since it 
is the product of the number of stars by the number of 
ionizing photons produced by one star. This degeneracy 
can be explored by changing the luminosity of the indi- 
vidual stars and by adjusting the corresponding value of 
the number of stars. 

If one changes the luminosity of the individual stars, 
the effective temperature of each star must be changed 
in order to recover the [Nil] 121.7/rni/[Nm] 57.3/im ratio. 
All the fluxes are then reproduced as in the adopted model 
presented in Table || within a few percent. Fig. || shows in 
an HR diagram the locations of the available CoStar mod- 
els (crosses) and the four models retained and discussed 
here (filled diamonds). The number of stars needed to re- 
produce the 2 cm flux density is given for each model. 
Although the range in effective temperature seems small 
(from 30 to 35 kK), it is large if one considers the strong 
constraint from the [Nil] 121.7/im/[Nni] 57.3/xm ratio (see 
Fig. |2|). Along the track between the four models with 1.5, 
3.3, 8.6 and 19 stars, the stellar age varies approximatively 
from 2.8, 3.5, 4.0 to 1.6xl0 6 yr. Whatever the exact num- 
ber of stars involved in the ionization of G29.96, we see 
that these stars occupy a rather small range in effective 
temperature (30 to 35 kK) and age (1.6 to 4 Myr). The 



obtained ages are quite old compared to "classical" expec- 
tations for UCHII regions. 

We cannot find a satisfying model corresponding to one 
single star, as the luminosity will then overstep the CoStar 
models grid and enter the post main sequence and/or Wolf 
Rayet (see Fig. ||). Once the temperature is derived from 
the diagnostic lines ratio [Nil] 121.7/xm/[Niii] 57.3/im, we 
used the most luminous star available and multiply its flux 
by a factor 1.5 to reproduce the radio flux ^. As we know 
from NIR observations (see Sect. ^^) that only one star is 
the primary source of ionization, we think our model with 
1.5 star is better. We can interpret the value of 1.5 star as 
a consequence of mixing one main ionizing star with one 
or more lower luminous stars. 



5.6. Dependence on atmosphere models 

The derived parameters of the ionizing source depend 
on the adopted atmosphere models. Given the limited 
amou nt o f constraints available on the ionizing fluxes (cf. 
Sect. 3.2) and the potential uncertainties of the CoStar 



models especially f or cool stars with weak winds ( Schaerei 
|fe de Koter 1997| ), we have also tested other non-LTE 
model atmo spheres. A full description is given in Morisset 
|ct al. (2002| ). Here we summarize the main effects. 

We have used the recent line blanketed models of 



Pauldrach et al. (200l[ ) and test calculations for O stars 
based on the c omoving frame code CMFGEN of Hillici 
|fc Miller (1998 ) which both include stellar winds. Spectra 
from the f ully blanketed plane par allel non-LTE TLUSTY 
models of Hubeny fe Lanz (1995 ) were also kindly made 
available to us by Thierry Lanz. The comparison of the 
predicted IR fine-structure line ratios with observations 
from the sample of Paper I & II shows that a consistent fit 
of all four ratios ([N m]/[N n], [Ar m]/[Ar n], [S iv]/[S ill], 
[Ne in]/[Ne n]) within a factor of two is only obtained 
with the CoStar models. The scatter in the effective tem- 
perature determined with the CoStar models is due to the 
[Ar ili]/[Ar il] ratio, which have a similar ionization poten- 
tial than [N III] (29.6 and 27.6 eV respectively), showing 
a potential problem in either the observed line fluxes, the 
attenuation correction process, or in the atomic data. The 
other three excitation ratio, tracing ionizing photons be- 
tween 27.6 and 40.1 eV, are in a very good agreement, as 
also seen for G29.96. 

Using the extreme excitation ratio of [N ni]/[N n] and 
[Ne m]/[Ne n], we estimate effective temperatures of ~ 
32 -35, 33-38 and 34-38 k K, using the models CMFGEN 
of [Hillicr fc Miller (1998"| ), TLUSTY of [Hubeny fc Lanz 
|(1995D , and WM- Basic of [Pauldrach et al. (200l| ) re- 



spectively, while the same ratio leads to lower and less 
scattered effective temperatures (29.5-30.5 kK) using the 
CoStar models (see fig. ||). 



In practice the SED we used in the model is a weighted 



sum of D4 ( 70 %) and D5 (30 %) models from Schaerer fc de 



Koter (1997), i.e. a 2.9 x 10 6 years old 07 star of 50 M G 
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5.7. Other constraints and comparison with earlier 
studies of the exciting star 

Various estimates of the properties of the dominant ioniz- 
ing star of G29.96-0.02 have been obtained from the fol- 
lowing I observations / methods: 



duce the observed A^yc for a temperature range between 
~ 30 and 46 kK, depending on the evolutionary state. 

3) Given the obvious importance of small number 
statistics for the number of massive stars observed in 



the cluster associated with G29.96 (sec Afflcrbach et al 



1997; Pratap et al. 1999) standard evolutionary synthe- 



1. the bolometric luminosity L, estimated from IR or 
multi-wavelength observations 

2. the photon flux in the Lyman continuum A^yc assum- 
ing a single star, 

3. the ratio N-^ yc /L for a stellar cluster, 

4. IR line ratios through nebular modeling, 



sis models cannot be used for comparisons of N-L yc /L 
(Cerviho et al. 200C). However, an analysis of the resolved 
stellar content provides further insight. As discussed by 



Afflcrbach et al. (1997 ); Pratap et al. (199E ) we have taken 



the objects with H—K > 1. as cluster members. Assu ming 
a mean extinction Ah — 3.6 ( Afflcrbach et al. 1997) and 



5. the|He + abundance derived from radio recombination 
lines, 

6. near-IR photometry, and 

7. K-band spectral classification of the central source. 

The estimated spectral types quoted in the literature reach 
from 03 to 09.5 (e.g., Lacy et al. 1982; Doherty et al 



using the synthetic photometry of Lcjcunc fc Schaerer 



(2001 ) we have determined from the H-band magnitude 



1994; Simpson et al. 1995; Afflcrbach et al. 1997; Watson 



et al. 1|99Y; Hapcr et al. 2002). However, these estimates 



rriH the luminosity of the individual stars assuming all 
members to be on the same isochrone with ages between 
~ and 4 Myr. From this we derive the fraction of L 
provided by the ionizing star, which is found to be ~ 70- 
50 % for ages 0-4 Myr. A somewhat smaller fraction (70- 
30 %) is found using mx (and Ak = 2.14). This quan- 
titative estim ate confirms the expectations of Afflerbacfj 



are partly based on incompatible hypothesis such as dit- 
ferent assumptions on the source distance and different 
spectral type vs. T c g calibrations. Furthermore only un- 
evolved zero age main sequence (ZAMS) stars were con- 
sidered in most cases, in conflict with recent evidence 



(Afflcrbach et al. 1997, this paper). For these reasons we 
here briefly rediscuss these estimates using also up-to-date 
stellar models and homogeneous assumptions. In particu- 
lar all observational constraints are scaled to a distance of 
6 kpc (cf. Sect. |2.3| ). As T G ff is the physically most impor- 
tant parameter - independently of the exact spectral-type 
vs. T c ff relation - we essentially derive the constraint on 
this parameter. 

1) The total bolometric luminosity of G29.96 obtained 
from the 12-100 (j,m IRAS flux and its overall SED in 
an arcminutc sized region is log L/Lq ~ 5.90 (Paper I, 
Afflcrbach et al. 1997). Likely the major fraction of it 
is due to the main ionizing source (cf. |Afflerbach et al" 



1997, and below). Using the calibrations of [gchmidt-Kaler 
(1982| ) yields the following T c g: ~ 44 kK (fo r luminosity 



class V ), 41 kK (LC III), 37 kK (LC I). From |Vacca et al 
p996| ) one obtains: ~ 48 kK (LC V), 45 kK (LCIII), 36 
kK (LC I). A very wide range of T cS 20 kK to < 50 
kK) is allowed for main sequence stars with the given lu- 
minosity, as shown in Fig. |^. These T c g represent upper 
limits as other stars contribute to the total bolometric (cf. 
below) . 

2) The ionizing photon flux derived for G29.96 from 
radio emission is log(A^L y c ) ~ 49.-49.14 s" 1 (Fey et al. 
1995| ; Kim & Koo 2001). A somewhat higher value of 



LycJ 



49.29 was derived by Afflerbach et al. (1997) 



log(JV] 

fro m the extinction cor rected Brackett-7 map. Based on 
the Vacca ct al. (1996 ) calibrations this corresponds to 
T eS ~ 40-43 kK (LC V), 35-38 kK (LC III), and 30-32 



kK (LC I). S imilarly, the stellar models of pchacrer fc de 
Kotcr (1997 ) (based on the tracks shown in Fig. |3|) repro 



et al. (1997) of a contribution of at least 50 % from the 
ionizing star to L. Correcting for ~ 50 % of L due other 
cluster members and assuming that one star dominates 
the ionization, we obtain a revised N-^ yc / L of the ionizing 
star which should be comparable to predictions for single 
stars. The comparison with the stellar models used earlier 
yields T c g between ~ 31 and 38 kK. 

4) |Simpson ct al. (1995| ) and [Afflerbach et al. (1997| ) 
used line measurements from KAO and photoionization 
models including plane parallel LTE Kurucz model at- 
mospheres. Their analysis (method 4) yields T c g = 35.7 
and 37.5 kK respectively, rather similar to our values de- 
rived with different atmosphere models, but larger than 
the value obtained with the CoStar atmosphere models. 
The main difference with their result likely originates from 
the use of more sophisticated model atmospheres. 

5) The He + /He ionization fractions derived for the 
best model presented in Table || are 50 % (35 %) for 
Component 1 (2) respectively, slightly lower than the val- 
ues obtained by [Kim fc Koo (200l[ ): 68 to 76 %. Using 
hottest stars with CMFGEN at 33 kK and WM-Basic at 
36 kK atmosphere models, we found He + /He to be 60 % 
(48%) and 77 % (68 %) resp ectively, in better a greement 
with the value obtained by Kim & Koo (2001) (but see 
the discussion on the X i+1 /X* ratios for Ar, Ne, N and S 
in Sect. j5~6[ ). 

6) From photometric observations and constraints on 
the total luminosity of G29.96 ( [Afflerbach et al. 1997j ) de- 
rive an allowed temperature range for the ionizing star of 
T eS ~ 28-37 (23-43) kK at for 1 (3) a uncertainties valid 
for source distances between approx. 5-9 kpc. Our above 
analysis of the cluster photometric data, taking the the 
contribution of all individual objects to L into account, 
yields consistency only for ages ~ 3-4 Myr. Despite this, 
the permitted T g range based the H or K band data re- 
mains fairly large, and essentially identical to the above 
values. 
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7) Watson & Hanson (1997) obtained the first K-band 
spectrum of the ionizing star of G29.96, whose spectral 
type was found between 05 and 08 (luminosity class un- 
determined; cf. Watson & Hanson 1997), based on the 
presence of He II absorption, and C iv and N in emission 
lines. They note, however, that a 07 or 08 spectral type 



van Burcn & Mac Low 1992; Lumsdcn & Hoare 1996) 



Champagne flow models (see e.g., Yorke et al. 1983), re- 
sulting from the expansion of an H II region into a molec- 
ular cloud exhibiting a density gradient, are also able to 
reproduce the cometary morphology. These models were 



applied more recently to G29.96 (Fey et al. 1995 



would require some enhancement of the C IV and N III & Hoare 1996, PMB99) and were found to give results 



Lumsdcn 



features - attributed to a higher metallicity - compared 
to "normal" objects of these types. While the recent VLT 



spectrum presented in the preliminary work of Kaper et al 



(2002[[ 



appears to be consistent with the data of Watson 



& Han json (1997| ), the former authors deduce a spectral 



type as early as 03 based on the presence of the C iv and 
N in emission lines. From this it appears that a m ore de- 
tailed analysis of the data of Kaper et al. (2002 ) should 
be awaited before more firm conclusions on the spectral 
type of G29.96 can be drawn. 

In any case, given the unknown luminosity class the fol- 
lowing temperature ranges are obtained for 05-08 (03): 
- 38.5-46 kK (51 kK) for LC V, intermediate val ues fo 
LC III, and ~ 36-45 kK (50 kK) for LC la using the |Vacca 



et al. ( 1996 ) compilation based on analysis using pure H- 
He atmosphere models. Recent fully line blanketed non- 
LTE calculations including stellar winds show, however, 
that - as already suspected earlier - the T c g scale of O 



stars must be cooler (e.g., Fullerton et al. 2000; Martins 
et al. 2|002|). 



The models of Martins et al. (2002 ) yield a reduction 
of T c ff by 4 to 1.5 kK for 03-09.5 dwarfs compared to 
the Vacca et al. (1996 ) scale, and larger reductions are 
expected for giants and supergiants. Taking these effects 
into account we estimate for 05-08 types T ff ~ 36-43 
kK for LC V and - 33-40 kK for LC la. 

Combining the available data it appears that the p re- 
liminary spectral classification by Kaper et al. (2002 ) is 
the only information which is incompatible with most 
other constraints (points 3-6, possibly also 1 and 2). Good 
consistency is obtained, however, from the intersection of 
the above constraints 1) to 6), yielding an allowed T c g be- 
tween ~ 31 and 37 kK, overlapping with the spectral type 
derived by Watson & Hanson (1997). We thus conclude 
that overall the parameters derived from our photoioniza- 
tion modeling are compatible with all the available obser- 
vational data. 



5.8. Implications of the advanced age of G29.96 

However, an age of ~ 3 x 10 6 years for the star is very high 
compared to the expected dynamical lifetime of UCHII re- 
gions (5 x 10 5 years, see e.g., Wood & Churchwell 1989a, 



based on the number of UCHII regions in the Galaxy and 
their expected lifetime). Two main models have been de- 
veloped to explain the cometary morphology which is com- 
mon for UCHII regions. Models of stellar- wind bow shocks 



(see e.g., |Mac Low et al. 199l| ), due to an O star mov- 
ing supersonically through a molecular cloud, were first 



more consistent with the radio observations. 

It is important to note that assuming a reasonable 
projected proper motion of 1 km/s, the star should have 
moved away from its birth place by about 3 pc (1.75 ar- 
cmin) in ~ 3 x 10 6 years. This rules out the Champagne 
flow model as a complete description of G29.96 and 
strongly favors the random meeting of an older star with 
an interstellar cloud. The ionizing star may also have left 
its birthplace, irradiating molecular gas further out which 
could still be part of the larger parental cloud from which 
it was formed. 

5.9. Reliability of the abundances determination 

The determination of the elemental abundances from the 
infrared fine-structure lines depends on many physical pa- 
rameters, such as the filling factor, which are poorly con- 
strained. Nevertheless, we can assert that there are two 
groups of elements. On one hand, oxygen and nitrogen, 
whose lines, all observed by the LWS spectrometer, are 
mostly emitted by the extended component 1, due to their 
low critical densities. Uncertainties in the attenuation cor- 
rection and then in the Bra line flux by, e.g., 10 % leads 
to an uncertainty on the N and O abundances of 25 to 



30 % (see Sect. |5jJ). 

The elements neon, argon and sulfur group, whose lines 
are observed by the SWS spectrometer (as the Hi lines) 
with all the subsequent aperture corrections, are emitted 
by both components. The pres ence of high density clumps 
(filling factor of 0.1 - see Sect. 5^2) in the core will lead to 
abundances two times higher than what we determined in 
the presented model for the Ne, Ar, S group. 

Whatever the uncertainties could be on the filling fac- 
tor, the geometry of the source, the attenuation or the 
actual value of the radio emission, the determination of 
the abundance ratios in each group are robust: the N/O 
ratio on one hand, and the Ne/Ar, Ne/S and Ar/S ratios 
on the other. 

Tabled compares the abundances determined here and 
by [Hcrtcr et al. (198l| ), |Simpson et al. (1995| ), [Afflcrbacb 



t al. (1997) and Paper II. The solar abundances from 
Grevesse fc Sauval (1998|) are also given. [Afflcrbach et al 



(1997) used thc |5impson et al. (1995 ) observations to mod 



studied and applied to G29.96 ( |Wood fc Churchwell 199l| ; 



eled G29.96, but with a core/halo description. They both 
made semi-empirical models (using icf 's) and adopted an 
higher effective temperature (~ 36 kK, see discussion in 
Sect. |5.7|) . The method used in Paper II is semi-empirical, 
based on the same observed line fluxes as the present work. 
For those previous works, we give the values effectively 
determined for G29.96 and the values obtained using the 
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Table 3. Abundances determined by previous authors and in the present work. 



Element 


Herter et al. (198]) 


Simpson et al. (1995) 


Afflcrbach et al. (1997) 


Paper II 


This work 


Solar^ 




G29.96 


G29.96 


4.5 kpc 1 


G29.96 


4.5 kpc 1 


G29.96 


4.5 kpc 1 






N/H (lO" 4 ) 




2.3 


1.8 


1.8 


1.2 


1.9 




2.0 


0.8 


O/H (lO" 4 ) 




8.5 


6.6 


5.6 


7.3 


5.1 




4.6 


6.8 


Ne/H (10~ 4 ) 


2.7 


2.6 


1.8 






2.5 


2.2 


1.7 


1.2 


S/H (HT 5 ) 


3.2 


1.9 


1.6 


2.2 


1.8 


0.8 




2.2 


2.1 


Ar/H (10~ 6 ) 


23. 










4.8 


5.0 


5.0 


2.5 


N/O 




0.27 


0.27 


0.32 


0.17 


0.37 


0.33 


0.43 


0.12 


Ne/S 


8.4 


13. 


11. 






36. 




7.5 


5.7 



Values obtained applying the gradients from the corresponding authors at the galactocentric distance of G29.96 



from Grevesse & Sauval (1998) 



abundance gradient law they found, applied at the posi- 
tion of G29.96: 4.5 kpc from the galactic center. 

The set of abundances, exepted for oxygen, shows that 
G29.96 is overabundant compared to the solar values, in 
agreement with its inner position in the Galaxy. The abun- 
dances determined in the present work are compatible 
with the previous determination within a factor of 2 , ex- 
cept with the Ar/H ratio from Herter et al. (1981 ) and 
S/H ratio from Paper II. 

The determination of the sulfur abundance relative to 
hydrogen in Paper II is very different from what the previ- 
ous authors and the present work found. From the results 
presented in Table || we see that the [S in] 18.7/zm and 
[S iv] 10.5/xm lines on which the sulfur abundance is based 
in Paper II are mostly emitted by the extended compo- 
nent 1. The effect of the finite aperture size of the SWS 
instrument is crucial in this case. As there was no correc- 
tion for this effect in Paper II, the sulfur emission and its 
abundance are underestimated. 



6. Conclusion 

This paper presents a detailed model of the compact 
Galactic Hn region G29.96 for which high quality data 
imagery and spectroscopy is available at both infrared 
and radio wavelengths, including recent ISO observations. 
The model, which is based on the photoionization code 
NEBU and state-of-the-art stellar atmosphere models, re- 
produces most of the observations, with the exception of 
a few points known to be less accurately measured. The 
radio and infrared data on G29.96 arc best reproduced 



effective temperature scale of O stars ( Martins et al. 2002 ) 
are taken into account. 

We explored the effect of varying the different model 
parameters on the predictions. The main sources of un- 
certainty in determining the abundances are the fluxes of 
the hydrogen recombination lines and the geometry of the 
dense compact core. 

The derived elemental abundances are in agreement 
with the lowest values found in previous studies. The most 
robust results are N/O and Ne/S which are 3.5 and 1.3 
times the solar values, respectively. 

From the reanalysis of the different available observa- 



tional constraints (see Sect. 5.7) it is not surprising that 
several earlier studies reached apparently conflicting re- 
sults on the spectral type or T e g of the ionizing source of 
G29.96. This is mostly due to the following facts. First, 
proper photoionization models must account for the de- 
pendence of nebular lines on several parameters, includ- 
ing the ionization parameter, geometry etc. Second, con- 
sistent predictions from stellar models regarding ionizing 
fluxes, T c g scales etc. should be used taking into account 
recent progress made with fully line blanketed non-LTE 
atmosphere models. Last, but not least, as the available 
K-band spectroscopy does not allow a determination of 
the luminosity class, allowance should be made for possi- 
ble variations of the evolutionary state of ionizing sources 
in compact Hn regions. In view of these issues one may 
question whether mid-IR analysis of compact H II regions 
yield intrinsically different estimates of T e g compared to 



spectral types or other constraints as suggested by Hanson 
et al. (2002|). Our detailed analysis of G29.96 indicates the 



by a 2-density component model nebula, with a diffuse 
(riff ~ 600 cm -3 ) extended (1 pc) halo surrounding a 
dense (n# ~ 50000 cm~ 3 ) compact (0.1 pc) core. 

Using CoStar stellar atmosphere models we derived an 
effective temperature of ~ 30^ kK. Adopting more recent 
non-LTE line blanketed atmospheres with stellar winds, a 
somewhat higher T e ff ~ 32-38 kK is found. This tempera- 
ture range is compatible with all observational constraints. 
For T e ff ~ 33-36 kK compatibility is also obtained with 



the K-band spec tral type 05-08 determined by ( |Watson 



fc Hanson 1997 ) when recent downward revisions of the 



opposite. 

The age of the ionizing star required by our model is 
3 x 10 6 yr, much older than the expected lifetime of UCHII 
regions. This could indicate that G29.96 is not excited by 
a bona fide young massive star. Instead the ionizing star 
creating today the Hil region G29.96 may have left its 
birthplace, exciting gas further out. This matter could still 
be part of the larger parental cloud from which the stellar 
cluster associated with G29.96 was formed. 

Acknowledgements. We thank the referee for useful questions 
and comments. CM thanks D. Pequignot for helpful discus- 



C. Morisset et al.: A photoionization model of G29.96-0.02 



13 



sions on photoionization models and Ryszard Szczerba for dis- 
cussions on dust. DS thanks Alan Watson, Margaret Hanson, 
and Yuri Izotov for useful discussions and Margaret Hanson 
for sharing data before publication. We thank Thierry Lanz 
for sending us line blanketed TL USTY model atmospheres be- 
fore publication, and John Hillier for making his atmosphere 
code CMFGEN available. 

References 

Aannestad, P. A. & Emery, R. J. 2001, A&A, 376, 1040 
Affierbach, A., Churchwell, E., Acord, J. M., Hofner, P., Kurtz, 

S., & Depree, C. G. 1996, ApJS, 106, 423 
Affierbach, A., Churchwell, E., Hofner, P., & Kurtz, S. 1994, 

ApJ, 437, 697 

Affierbach, A., Churchwell, E., & Werner, M. W. 1997, ApJ, 
478, 190 

Cervino, M., Luridiana, V., & Castander, F. J. 2000, A&A, 
360, L5 

Claussen, M. J. & Hofner, P. 1995, in Bulletin of the American 

Astronomical Society, Vol. 27, 1440 
Doherty, R. M., Puxley, P., Doyon, R., & Brand, P. W. J. L. 

1994, MNRAS, 266, 497 
Draine, B. T. & Laor, A. 1993, ApJ, 402, 441 
Fey, A. L., Gaume, R. A., Claussen, M. J., & Vrba, F. J. 1995, 

ApJ, 453, 308 

Fullerton, A. W., Crowther, P. A., De Marco, O., Hutchings, 
J. B., Bianchi, L., Brownsberger, K. R., Massa, D. L., 
Morton, D. C, Rachford, B. L., Snow, T. P., Sonneborn, 
G., Tumlinson, J., & Willis, A. J. 2000, ApJL, 538, L43 

Garay, G., Rodriguez, L. F., Moran, J. M., & Churchwell, E. 
1993, ApJ, 418, 368 

Garcia-Lario. 1999, 2nd Post-Operations Cross-Calibration 
Status Report, Tech. Rep. SAI/1999-083/Rp, ISO Data 
Center, http:/ /www.iso.vilspa.esa.es/users/expLlib/ISO/ 

Giveon, U., Sternberg, A., Lutz, D., Feuchtgruber, H. , & 



Martin s, F., Schacrer , D. 



see istro-ph/0111233 



Hillier, D. 2002, A&A, in press, 



Pauldrach, A. 2002, ApJ, in press, see a,stro-ph/01 10366 



Grevesse, N. & Sauval, A. J. 1998, Space Science Reviews, 85, 
161 

Hanson, M., Luhman, K., & Rieke, G. 2002, ApJS, in press 
Herter, T., Heifer, H. L., Forrest, W. J., McCarthy, J., Houck, 

J. R., Willner, S. P., Puetter, R. C, Rudy, R. J., Soifer, 

B. T., & Pipher, J. L. 1981, ApJ, 250, 186 
Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 
Hubeny, I. & Lanz, T. 1995, ApJ, 439, 875 
Kaper, L., Bik, A., Hanson, M., & Comeron, F. 2002, in ESO 

workshop on The Origins of Stars and Planets: the VLT 

view, ed. J. Alves & M. McCaughrean (ESO Astr. Symp.), 

in press 

Kennicutt, R. C, Bresolin, F., French, H., & Martin, P. 2000, 

ApJ, 537, 589 
Kim, K. & Koo, B. 2001, ApJ, 549, 979 

Kurtz, S. E., Watson, A. M., Hofner, P., & Otte, B. 1999, ApJ, 
514, 232 

Lacy, J. H., Beck, S. C, & Geballe, T. R. 1982, ApJ, 255, 510 
Lejeune, T. & Schaerer, D. 2001, A&A, 366, 538 
Lumsden, S. L. & Hoare, M. G. 1996, ApJ, 464, 272 
Mac Low, M., van Buren, D., Wood, D. O. S., & Churchwell, 

E. 1991, ApJ, 369, 395 
Martin-Hernandez, N., Peeters, E., E. Damour, F., Cox, P., 

Roelfsema, P., Baluteau, J.-P., Tielens, A., Churchwell, E., 

Jones, A., Kessler, M., Mathis, J., Morisset, C, & Schaerer, 

D. 2002, A&A, 381, 606, Paper II 



Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & 

Charbonnel, C. 1994, A&AS, 103, 97 
Morisset, C. & Pequignot, D. 1996, A&A, 312, 135 
Morisset, C, Schaerer, D., & Bourret, J.-C. 2002, A&A, in 

preparation 

Oey, M. S., Dopita, M. A., Shields, J. C, & Smith, R. C. 2000, 
ApJS, 128, 511 

Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 2001, 
A&A, 375, 161 

Peeters, E., Martm-Hernandez, N., E. Damour, F., Cox, P., 
Roelfsema, P., Baluteau, J.-P., Tielens, A., Churchwell, E., 
Jones, A., Kessler, M., Mathis, J., Morisset, C, & Schaerer, 
D. 2002, A&A, 381, 571, Paper I 

Pequignot, D., Ferland, G., Netzer, H., Kallman, T., 
Ballantyne, D., Dumont, A.-M., Ercolano, B., Harrington, 
P., Kraemer, S., Morisset, C, Nayakshim, S., Rubin, R., 
& Sutherland, R. 2002, in Photoionized Plasma, Lexington 
2000, ed. G. Ferland, in press 

Pratap, P., Megeath, S. T., & Bergin, E. A. 1999, ApJ, 517, 
799, pMB99 

Schaerer, D. 2000, in ASP Conf., Vol. 221, Stars, Gas and Dust 
in Galaxies: Exploring the Links, ed. D. Alloin, G. Galaz, & 
K. Olsen, 99 

Schaerer, D. & de Koter, A. 1997, A&A, 322, 598, 



http : / / webast . ast.obs-mip.fr /people / schaerer / 



Schaerer, D. & Stasiriska, G. 1999, in The Universe as seen 

by ISO, ed. P. Cox & M. Kessler, ESA Special Publications 

series (SP-427), 751 
Schmidt-Kaler, T. 1982, in Landolt-Bornstein, New Series, 

Group VI, ed. K. Schaifers & H. Voigt, Vol. 2b (Berlin: 

Springer) 

Sellmaier, F. H., Yamamoto, T., Pauldrach, A. W. A., & 

Rubin, R. H. 1996, A&A, 305, L37 
Simpson, J. P., Colgan, S. W. J., Rubin, R. H., Erickson, E. F., 

& Haas, M. R. 1995, ApJ, 444, 721 
Stasiriska, G. & Schaerer, D. 1997, A&A, 322, 615 
Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 

914 

van Buren, D. & Mac Low, M. 1992, ApJ, 394, 534 

Watarai, H., Matsuhara, H., Takahashi, H., & Matsumoto, T. 

1998, ApJ, 507, 263 
Watson, A. M., Coil, A. L., Shepherd, D. S., Hofner, P., & 

Churchwell, E. 1997, ApJ, 487, 818 
Watson, A. M. & Hanson, M. M. 1997, ApJL, 490, L165 
Wood, D. O. S. & Churchwell, E. 1989a, ApJ, 340, 265 
— . 1989b, ApJS, 69, 831 
— . 1991, ApJ, 372, 199 

Yorke, H. W., Tenorio-Tagle, G., & Bodenheimer, P. 1983, 
A&A, 127, 313 



